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Reaction-in-Flight Neutrons as a Signature for Shell Mixing in NIF capsules 
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We present analytic calculations and results from computational simulations showing that 
reaction-in-flight (RIF) neutrons act as a robust indicator for mixing of the ablator shell mate- 
rial into the fuel in DT capsules designed for the National Ignition Facility. The sensitivity of RIF 
neutrons to hydrodynamical mixing arises through the dependence of RIF production on charged- 
particle stopping lengths in the mixture of DT fuel and ablator material. Since the stopping power 
in the plasma is a sensitive function of the electron temperature and density, it is also sensitive to 
mix. RIF production scales approximately inversely with the degree of mixing taking place, and the 
ratio of RIF to downscattered neutrons provides a measure of the mix fraction and/or the mixing 
length. For sufficiently high-yield capsules, where spatially resolved RIF images may be possible, 
neutron imaging could be used to map RIF images into detailed mix images. 
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At the National Ignition Facility, reaction-in-flight 
(RIF) neutrons are expected to make up to about 1.5% 
of the total neutron production in high-yield capsule im- 
plosions. In this paper, we show how RIF production 
is sensitive to mixing of the capsule shell material into 
the DT fuel and how this sensitivity can be used as a 
probe of such mixing. RIF neutron production is af- 
fected by mix through a direct dependence on charged- 
particle stopping. Because mixing causing changes in 
electron temperature and density, it directly affects the 
charged-particle range in the plasma. Energetic ions un- 
dergo fewer reactions-in-flight if they are stopped more 
quickly, so the number of RIF neutrons decreases with a 
decrease in the charged-particle range. As discussed be- 
low, our radiation-hydrodynamic-burn simulations con- 
firm the expectation that the magnitude of the RIF neu- 
tron spectrum (relative to the total neutron production) 
scales directly with mix. This is in contrast to the other 
components of the neutron spectrum, which, to first or- 
der, simply scale up and down with the yield. 

RIF neutrons are produced by a two-step process. In 
the first step, a primary 14.1 MeV DT neutron knocks a 
triton or deuteron up to a spectrum of energies from zero 
to more than 10 MeV for tritons. In the second step, the 
energetic knocked-on ion undergoes a DT reaction with a 
thermal ion, producing a neutron above the primary 14.1 
MeV peak. The continuous spectrum of neutrons pro- 
duced in this way constitute the RIF neutron spectrum. 
The spectrum of knock-on particles after transport away 
from the point of production has a similar shape, but is 
distorted by energy loss in the DT fuel [J, 0] • 

A second mechanism for production of RIF neutrons 
begins with energetic a particles born with 3.5 MeV of 
kinetic energy in the reaction d+t — > n + a. These a par- 
ticles collide with ions, producing another distinct com- 
ponent of high-energy knock-ons. RIF neutrons produced 
by reactions of these high-energy particles contribute to 
the RIF spectrum up to about 17.6 MeV. Their produc- 
tion involves two powers of the charge-particle stopping 
power, one for the alpha particle and one for the knocked- 
on triton or deuteron. For this reason, a-induced RIF 
production shows a different sensitivity to mix than nor- 
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FIG. 1: Differential cross section for knock-on production by 
elastic scattering of 14.1 MeV neutrons on deuterons and tri- 
tons, as a function of the energy of the scattered ion. These 
spectra represent the energy distribution for knock-on ions at 
the point of production. 



mal RIFs. Because the charged particle stopping is short 
in the DT fuel of a NIF capsule, the ratio of alpha- 
induced to neutron-induced knock-on ions is small. For 
the rest of this paper, we concentrate on neutron-induced 
RIFs only. 

The mix dependence of RIF production comes into 
play in the transport and energy loss of the energetic 
deuteron or triton, after it is knocked on. This charged- 
particle transport is sensitive to the plasma conditions [H, 
2]; in the relevant temperature regime (keV), charged- 
particle stopping power is dominated by scattering of ions 
with electrons in the plasma and behaves like0, 0] 

dE n„ , . , i in 

where 9, n e are the electron temperature and number 
density and E is the kinetic energy of the ion; log(A) 
is the Coulomb logarithm that describes the interplay 
between the long- and short-range Coulomb interac- 
tions, and c is a constant very nearly independent of 
plasma conditions. In principle, the energy distribution 
of charged-particles at a point distant from their point 
of production can be calculated directly in terms of the 
distribution at the point of production, using this stop- 
ping power. The rate of RIF reactions at any point is the 
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product of the deuteron/triton density at that point, the 
cross-section for d+t —> a+n, and the flux of transported 
knock-on ions at that point. The total RIF production 
over the whole capsule is the integral of this local rate and 
clearly depends directly on the charged-particle stopping 
length and the processing of the knock-on production 
spectrum by charged-particle transport. 

Formally, the local RIF production rate is defined by 
an integral over charged-particle paths, from the knock- 
on production point to the RIF reaction point. Let y 
be the point of production for a knock-on ion, and let 
x be the RIF reaction point. The path is the segment 
connecting y to x. Suppressing the energy dependence, 
we can write 
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where </>„ is the primary neutron flux and tpko is the re- 
sulting knock-on fluence (per unit initial particle), and 
CTfco and <Jdt are cross sections determining the knock-on 
and d + t reaction rates. Evaluating eq. (fTJ) requires, in 
general, detailed simulations that include treatment of 
the charged-particle transport. 

Before turning to our simulations, we note that con- 
siderable insight into the physics of RIF production can 
be obtained by expressing the charged-particle transport 
in terms of average paths for neutrons and for knock-on 
particles. Without loss of generality, the total number of 
charged-particle knock-ons produced in the capsule can 
be written 

— ^primary {^dt O~ko)nt 

where -/V pr j mary is the number of 14.1 MeV neutrons from 
the primary d + t reactions; the notation (■ • •}„ is short- 
hand for the integral over neutron paths, averaged over 
all paths taken by the d + t neutrons produced in the 
capsule. Because the capsule is thin to neutrons, the 
average length of neutron paths is, up to dimensionless 
geometric factors, equal to the size of the system. In a 
system thick to neutrons, the average length would be 
equivalent to the neutron mean free path. Similarly, the 
number of RIFs is 
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where fidt = \ndt is the density of either D or T ions, 
depending on whether the knock-on particles are T or D 
ions, respectively. We assume equal populations of D and 
T for our capsules of interest, which leads to the explicit 
factor of \ in this definition. Clearly the lengths of knock- 
on particle paths scale directly with the charged-particle 
stopping length, since the stopping length is small com- 
pared to the size of the system, for NIF capsules. 

Continuing with the qualitative discussion, we can fac- 
tor the above expressions into products of averages, each 
with the correct dimensions, and each scaling appropri- 
ately with changes in the underlying plasma properties. 



TABLE I: Definition of the mix fraction / in terms of the ion 
and/or electron densities in the gas. 
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The number of RIF reactions taking place is then 

N-rif « -N w i ma , ry (Tko<Tdt n 2 dt R Iko (2) 
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2 n e clog(A) 

where R is the size of the system. In a limit where the 
plasma is thin to charged-particles, Nrif would scale 
with (pr) 2 or, equivalently, with n 2 dt R 2 . However, NIF 
capsule implosions of interest are not in this regime. 
As indicated above, when the charged-particle stopping 
length is short compared to the size of the system, Nrif 
scales with only one power of ndt- The second power of 
ndt only appears as a ratio to the electron density, and 
charge neutrality requires that the ratio of ions to elec- 
trons rii/n e is a constant dependent only on the effective 
charge of the ions; to first order, variations in this ratio 
only occur with changes in material composition, such as 
those which occur due to mixing. 

Therefore, we next consider the effect of mixing on 
RIF production in this simple formulation. Using the 
definitions of mix fraction / in Table |H Eqn. (J3]) implies 
that the change in RIF production, relative to a clean 
no-mix (/ = 0) situation, is simply 
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Qmix j g e i ec t r on temperature in the presence of mix 
and 8® T is the temperature in the unmixed fuel. Z e g 
is the number of free electrons introduced into the gas 
per mix atom. Going beyond these assumptions, which 
lead to a single representative charged-particle stopping 
length, requires realistic simulations, such as those pre- 
sented below. We note that the function describing the 
dependence of the RIF fraction on the mix fraction / is 
sometimes referred to as the transfer function. 

To further fix the parameters of this qualitative dis- 
cussion, we can consider models which couple the change 
in electron temperature to the mix fraction, thereby fix- 
ing the ratio in Eqn. ([¥]) explicitly in terms of the mix 
fraction /. The simplest model of this type, which is not 
far removed from the true physical situation, stipulates 
an adiabatic mixing process. In the adiabatic limit, in 
which heat is conserved, the mix-induced change in the 
electron temperature can be related to the electron tem- 
peratures in the original un-mixed shell and un-mixed 
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DT fuel regions by 

Qmix = (1 + Z eS )f0f M + 2(1 - f)9° T 

(1 + 3*)/ + 2(1-/) ' {> 

This change in temperature as a function of mix fraction 
is shown in Fig. [2j 
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FIG. 2: (color online) The temperature in the mix region 
as a function of mix fraction f. In these calculations, Z e g 
was taken to be 6. The unmixed shell material temperature 
was assumed to 300 eV, which is a typical value obtained in 
simulations. The different curves are for three different DT 
electron temperatures, as indicated in the legend. 

Combining eqs. (4) and (5), the expected shape of the 
RIF production rate with increasing mix fraction / is as 
shown in Fig. [3J As can be seen, the rate of suppression 
of the RIF production with increasing mix fraction is 
controlled by Z c ff. Of course, in reality it is unlikely that 
there would be a single mix fraction or 9™ zx throughout 
the capsule; rather, / = f(x,y,z) in cylindrical coordi- 
nates. For high-yield shots, the RIF production is ex- 
pected to be large enough to allow for spatially resolved 
RIF neutron images l.'i . Using Eqs. 4 and 5, (or, equiva- 
lently, Fig. 2), a RIF neutron image could be translated 
into a spatially resolved map of the mix fraction. In 
other words, the so-called transfer function determines 
the mapping needed to convert a two-dimension neutron 
image of RIF/Total(a:, y) into an f(x, y) mix image. Such 
a capability would open the possibility of detailed hy- 
drodynamical mixing experiments at the NIF, assuming 
that yields sufficient for a Uses of Ignition program are 
obtained. 

Equations [4] and [5] were derived under the assump- 
tions that complications from charged-particle transport 
and capsule dynamics could be ignored. To examine 
the changes in RIF production with increasing mix frac- 
tion under more realistic conditions, we ran a series 
of one-dimensional radiation- hydrodynamic-burn simula- 
tions. Hydrodynamic mixing was treated using the Scan- 
napieco and Cheng mix model Q, which is essentially a 
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FIG. 3: (color online) RIF production, compared to a clean 
calculation, as a function of the mix fraction / for different 
values of Z c g, the number of free electrons introduced into 
the gas per mix atom. The curves shown are derived from 
eqs. (TJ][5j) which makes the over-simplified assumption that 
the material properties do not vary in the RIF production 
region. The temperatures 9® T are the assumed temperatures 
of the electrons in the DT fuel without mix (f m ix = 0) . 



fluid interpenetration model. In this model, the degree of 
mixing is controlled by a single parameter a, which con- 
trols the dynamic mixing length. Analyses 6] of Omega 
experiments typically require a value of a ~ 0.06 in or- 
der to reproduce observed yields. In the present calcu- 
lations, we allowed mixing across both the shell/ice and 
ice/gas interfaces. We varied a between 0-0.16, but al- 
ways used the same value of a at the two interfaces. 
Charged-particle transport was calculated using an ef- 
fective diffusion model. We have also conducted tests 
of the basic physics of the RIF production mechanism 
outlined in this paper, under more simplified physical 
circumstances, but using highly accurate Monte Carlo 
charged-particle transport. More accurate Monte Carlo 
transport methods could change the results of the full- 
physics simulations quantitatively, but will not change 
the conclusions of the paper. 

Since the RIF production depends on the areal density 
of the fuel, the RIF dependence on mix is most simply dis- 
played by dividing out the areal density. For this reason 
we examined the ratio of the RIF neutrons to downscat- 
tered neutrons, since the latter are a direct measure of the 
fuel (pr) 0, [1] . Our definition of RIFs included all neu- 
trons above 15 MeV, while the downscattered neutrons 
included all neutrons in the energy range 10-12 MeV. The 
energy restriction on the latter range avoids complica- 
tions from neutrons produced in t+t and other reactions. 
Choosing different energy ranges for either the RIF or 
downscattered neutrons did not significantly change the 
sensitivity of RIFs to mix. The results are shown in Fig. 
HI where a steady drop is seen in the RIFs to downscat- 
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tered ratio as the mix is increased. The decrease in the 
RIF ratio scales approximately linearly with the mixing 
length, until the degree of mixing is sufficiently large to 
induce failure. 
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In summary, both analytic arguments and computa- 
tional simulations show that RIF neutrons act as a ro- 
bust signature for mix in NIF capsules. The RIF neu- 
trons could be measured either by neutron time-of-flight 
@ or by using radiochemical techniques [ToL Ull. [l2| . Neu- 
tron imaging could also be used to measure RIF in high- 
yield capsules, and such images could be used to ob- 
tain spatially resolved mix images. Finally, we note 
that the shape of the RIF spectrum is determined in 
a relatively straightforward way from the shape of the 
knock-on spectrum, and the subsequent energy loss of 
the knock-on charged-particles in the plasma. If a-RIFs 
are included, the lower energy portion of the RIF spec- 
trum (E n ~ 15 — 18 MeV) has a more complicated shape, 
and this can be used to extract additional information on 
mix and electron temperature. 



FIG. 4: The panel on the left shows the RIF to downscattered 
neutron ratio with increasing mixing of the shell material into 
the fuel. The panel on the right is the RIF to total neuron 
ratio, a is the parameter in the Scannapieco and Cheng model 
that controls the mixing length. Here the RIF neutrons are 
all neutrons above 15 MeV, the downscattered are neutrons 
in the energy range E n = 10-12 MeV, the total are all those in 
the range 0-30 MeV. 



[1] P.B. Radha, S. Skupsky, R.D. Petrasso, J.M. Sources, 

Phys. Plasmas 7, 1532 (2000). 
[2] D.G. Hicks, C.K. Li, R.D. Petrasso, F.H. Seguin, B.E. 

Burke, J. P. Knauer, S. Cremer, R.L. Kremens, M.D. 

Cable,, and T.W. Phillips, Rev. Sci. Instrum. 68, 589 

(1997). 

[3] C. K. Li and R. D. Petrasso, Phys. Rev. Lett. 70, 3059 
(1993). 

[4] L. S. Brown, D. L. Preston, and R. L. Singleton, Phys. 

Rep. 410, 237 (2005). 
[5] A.J. Scannapieco and B.L. Cheng, Phys. Letts. A 299, 

49 (2002). 

[6] D.C Wilson, A.J. Scannapieco, C.W. Cranfill, M.R. 

Clover, N.M. Hoffman, Phy.s Plasmas, 10, 4427 (2003). 
[7] D.C. Wilson, W.C. Mead, L. Disdie, M. Houry, J.-L. 

Bourgade, T.J. Murphy, Nucl. Inst. Meth. A 488, 400 



(2002). 

[8] H. Azechi, M.D Cable, R.O. Stapf, Laser Part. Beams, 
9, 119 (1991). 

[9] Glebov, V. Y. et al., Rev. Sci. Instrum. 77,10E715-l-7 
(2006). 

[10] A.C. Hayes, G. Jungman, J.C. Solem, P.A. Bradley, R. 

S. Rundberg, Mod. Phys. Lett. A 21, 1029 (2006). 
[11] CP. Grim, et al., Rev. Sci. Instrum. 79, 10E503 (2008). 
[12] P.A. Bradley, T.A. Bredeweg, S. Cowell, M.M. Fowler, 

G.P. Grim, G.M. Hale, A.C. Hayes, G. Jungman, A.L. 

Keksis, R.C. Korzekwa, A. Obst, R.S. Rundberg, P.G. 

Sanchez, J.C. Solem, J.B. Wilhelmy, Los Alamos Internal 

Report, LA-UR-07-5428, and in preparation. 
[13] G.P. Grim, in preparation. 



